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We consider a Markovian jumping process which is defined in terms of the jump-size distribution 
and the waiting-time distribution with a position-dependent frequency, in the difi'usion limit. We 
I^ ' assume the power-law form for the frequency. For small steps, we derive the Fokker-Planck equation 

and show the presence of the normal diffusion, subdiffusion and superdiffusion. For the Levy distri- 
bution of the step-size, we construct a fractional equation, which possesses a variable coefficient, and 
^^ solve it in the diffusion limit. Then we calculate fractional moments and define fractional diffusion 

'^H , coefficient as a natural extension to the cases with the divergent variance. We also solve the master 

QQ . equation numerically and demonstrate that there are deviations from the Levy stable distribution 

^vj ' for large wave numbers. 



PACS numbers: 02.50.Ey, 05.40.Fb, 05.60.-k 



I. INTRODUCTION 



The diffusion is called anomalous if the mean squared displacement of the Brownian particle does not rise with 
time linearly, as it is the case for the usual Brownian motion, but slower (subdiffusion) or faster (superdiffusion). 
. 1 There are many examples of physical systems which exhibit such anomalous transport; they involve complex systems, 
trt I disordered systems, semiconductors, polymers, glasses, turbulent plasma, and many others (for a review see [EHl)- 
^ , Physical situations in which the deviation of the linear dependence is expected comprise: a porous and inhomogeneous 
environment with traps and barriers, a coupling to a fractal heat bath via a random matrix interaction |3|], long-range 
'^ , and/or long-time correlations. Also the transport phenomena in dynamical hamiltonian systems often exhibit anoma- 
Ph ■ lous behaviour because the corresponding phase space can contain some regular structures which act as dynamical 
O , traps |^|5|. A trajectory sticks to such self-similar hierarchy of islands in chaotic environment and abide on regular 
paths for a long time. The emergence of the anomalous diffusion means that the traditional approach, involving 
the standard Fokker-Planck equation (FPE) with the diffusion constant, is no longer valid. Usually, the anomalous 
diffusion is attributed to memory effects, which are especially pronounced in non-Markovian descriptions, e.g. by 
K*' means of the generalized Langevin equation |g and the continuous-time random walk models (CTRW) 1] . In the 
^^ diffusive limit, those approaches resolve themselves to the fractional equations formalism 0, Qilg which contain the 
Riemann-Liouville operator as the fractional time-derivative; the step-size distribution is the Gaussian. As a direct 
consequence of the nonlocality in time and, therefore, of the lack of Markovian property, the uncoupled CTRW pre- 
dicts subdiffusive solutions in this case j^, llfl • However, the anomalous behaviour is possible also for the Markovian 
\^ I processes; they are described by FPE with the position-dependent diffusion term and physically can correspond e.g. 
^^ to the diffusion on fractal objects [ll| or to the Langevin equation with the multiplicative noise |l3|- In this paper 
we discuss another example of such Markovian process. 

A more general approach considers the Levy distribution which is defined, in its symmetric version, via the char- 
acteristic function 4>{k) = exp(— a/c**), where 0</x<2. li n ^ 2 the second moment is infinite and the generalized 
central limit theorem must be applied. Physically, Levy processes can refiect self-similarity: they are the stable solu- 
tions of the renormalization group method and are invariant under the scaling of position and time. In contrast to the 
Q ' Gaussian distribution, they include large fluctuations. The Levy distributions are present in many problems connected 
CJ . with various branches of since, including not only physics but also biology, economy, flnancial market research, etc. 
L| ' The diffusion equation for the Levy process involves a fractional derivative over the process value [a, [ij, ll4L [la | and 
• ^H , for /x = 2 it resolves itself to the ordinary FPE. Since the variance diverges for /i < 2 and the traditional description 
/\^ • of the diffusion process is no longer valid, one has to introduce new concepts, e.g. to study fractional moments [ig 
M ' or to restrict the integration in averaging to a finite box which scales with time in a prescribed way |l7l llq | . 
- - I Problems related to generalized diffusion equations, which contain either anomalous behaviour of the variance or 

infinite fluctuations, are the subject of the present paper. We deal with a jumping process which is Markovian, 
defined in terms of a jump-size distribution Q{x) and the waiting time distribution Pp{t). A peculiarity of the 
distribution Pp{t) consists in a position-dependent jumping rate. The process is defined in Sec. II. In Sec. Ill we derive 
the FPE as a small step limit of the master equation. Sec. IV is devoted to the fractional diffusion equation which 
is an approximation of the master equation in the diffusion limit for Q given by the Levy distribution; we solve this 
equation, discuss its diffusion properties and compare with numerical solution of the master equation. The results of 
our analysis are discussed in Sec.V. 



II. DEFINITION OF THE PROCESS 

The jumping process we are to deal with in this paper is a stepwise constant stochastic process x(t) defined in terms 
of two probabiUty distributions [l9| . The waiting time density distribution determines the time intervals r^ between 
consecutive jumps and it is assumed in the Poissonian form: 

Pp{t) = ^.(a;)e-''(-)^ (1) 

where the jumping rate r^(x) depends on the process value (the position). The size of the jumps is determined by a 
given normalized distribution Q(r = x — x'). Then the infinitesimal transition probability can be easily constructed 
and the master equation derived. It is of the form 

—p{x,t)^-v{x)p{x,t)+ I Q{x ~ x')v{x')p(x' ,t)dx . (2) 

Because of the variable jumping rate v{x)^ the above process can be regarded as a version the kangaroo process 
|2(]ll21| . The difference consists in the jump-size dependence of Q - in the kangaroo process Q depends only on the 
current position. Due to that property, the Eq.j^l) can describe transport phenomena [l9j. On the other hand, the 
process constitutes a special case of the coupled CTRW. Taking into account the position dependence of jumping 
rate is an important generalization of traditional random walk approaches. Such dependence is expected in many 
phenomena in which inhomogeneity of the environment cannot be neglected 22j. 

The stationary solution of Eq.(|3l can be easily obtained: P{x) = l/v{x). The normalization condition imposes 
restrictions on the function v^x): it must rise sufficiently fast in the infinity and sufficiently slow near zero. If, in 
turn, llv[x) has poles at some x, the stationary solution also exists and it is in the form of a combination of the delta 
functions. In the other cases P{x) does not exist. A special version of the stationary process, defined on the circle, 
exhibits long-time correlations and can serve as a model of the 1// noise |2^ . 

The general, time-dependent solution of the Eq.(|21) can be obtained by using the Laplace transforms. The formal 
expression for the Laplace transform of p(a;, t) is the following pL9| 

Po{x) 1 f v{xq)pq{xo)Q{x - Xq) 
Pix, s) = — — -— + — — --- / ■ — dxo+ 

1 .^ f iy{xa)po{xo) tt i^(a;i_i)(3(a;i_i - a;i_2) , 

^ ^ n — tth^) — '^-^ 



+ — — r^Z^ / : — 7 — ^Q(x~xk-i) 

s + t^(.T)^7 s + iy{xo) 



dxo, 



where po i^) stands for the initial condition. 

By multiplying the Eq.Q by x'^ and by integrating over x, one can obtain the equation which governs the time 
evolution of the variance. Assuming that Q has the vanishing mean and finite second moment, we yield the following 
result: 



d{x^)p 
dt 



^')qW)p- (4) 



The simple case vix) =const can be solved completely and the Eq.Q predicts the normal diffusion. However, if {x)q 
does not vanish, the diffusion becomes ballistic [l9j. For Q with divergent second moment, in turn, {x'^)p is infinite. 

III. SMALL JUMPS: THE FOKKER-PLANCK EQUATION 

The expression Q is difficult to handle and one has to resort to approximations. In the limit of small jumps, the 
process becomes continuous both in space and time and FPE may be a candidate for such approximation. In order 
to construct it, we apply the Kramers-Moyal expansion |24j of the master equation. In that method, one changes the 
integration variable in the Eq.Q by introducing the step size r = x ~ x' and, after the expansion of the function 
p{x — r, t)v{x — r) around r = x, one obtains the master equation in a form of the following series 

ap(x,t)_^(-i)v^„, ^d 



r- 



dt ^—' n\ dx 

n—l 



)q(— )"[p(x,t)-Kx)] (5) 



which is still exact. The approximation consists in neglecting all terms of the order higher than two. Obviously, all 
moments of Q must be finite. Finally, we obtain the FPE: 



dp{x,t) a^ d^[v{x)p{x,t)] 



(6) 



dt 2 9x2 

where u^ = {r'^)Q- Therefore, the jumping rate may be interpreted as the position-dependent diffusion coefficient 
The approximation is valid if the jumps are small and the function v{x) is smooth |25| : 



Q{r) « for r > (5 

i'{x + Ax) « v{x) for Aa; < 5, 



(7) 



where (5 is a small parameter. 

For ^{x) =const, Eq.lO takes the form of the ordinary diffusion equation which describes the Wiener process and 
is characterized by the normal diffusion. An interesting case is the power-law form of the jumping frequency: 



v{x) = 7|x|" 



(8) 



where is a constant parameter and 7 ensures proper units; in the foUowingwe assume 7=1. The diffusion coefficient 
in the form JSjl has been used to describe the diffusion on fractal objects [ll|, the transport of fast electrons in a hot 
plasma 26] , turbulent two-particle diffusion |27| . The FPE with that diffusion coefficient has been analyzed from the 
point of view of the first passage time in Ref. 28] . The FPE ® with v given by Eq. Q can be solved exactly |23 : 



p{x, t) = Ce 



\x\^ exp(- 



2\x[ 



7'^(2+eyH 



{aH/2)^ 



(9) 



where Ce 






0\ 2+" . Some values of 9 must be excluded. Since for 6 < —3 the Eq.Q has the stationary 

solution (in the form of the delta function) , the approximation by the FPE is not valid; the diffusion in this case is 
local. Moreover, for 6 G (—2,-1) the distribution (|5J| is not normalized. Therefore, we have either 6 e (—3,-2) or 
9>-l. 




FIG. 1: Time evolution of the FPE solutions p{x,t) calculated from the Eq.Q (solid line) and solutions of the master equation 
@ for the Gaussian form of the distribution Q (dots). Parameters are: 6 — 0.5 and a — 0.1. 



Fig.l presents the comparison of the FPE solutions p{x,t), calculated according to the Eq.®, with the solutions 
of the master equation J^Jl, obtained from Monte Carlo simulations of the stochastic trajectories. The agreement is 
good already for a = 0.1. 

The diffusion properties of the system follow directly from the solution ^. The mean squared displacement for 
> — 1 is given by the power-law function of time by the following formula 



(x'it)) 



Y(3±l) 
^ \2+e) 

^ \2+e) 



'-(2 
2 ^ 



(10) 



and the diffusion coefRcient is: 

p = Ml ^ i-«/(M+e) (i^oo). (11) 

Then we can distinguish three cases. For 9 G (—1, 0), V is infinite and the superdifFusion emerges. For > 0, P = 
and we get the subdifFusion. Finally, the normal diffusion takes place for 9 = Q. Therefore, the jumping process 
involves all kinds of diffusion. 

IV. THE FRACTIONAL DIFFUSION EQUATION 

Let us now assume the distribution Q in the form of the symmetric Levy distribution, defined by its Fourier 
transform: 

Q(A;)=exp(-X^|fc|^) [K > Q), (12) 

where 1 < /i < 2. In contrast to the Gaussian distribution, Q{x), corresponding to the Ea. (|12|l . has algebraic tails, 
Q{x) ~ |a;|~^~^ (|a;| — > oo), which makes long jumps very probable. In the diffusion limit fc — > 0, the Ea. (|12|l is given 
by Q{k) « 1 — K^\k\^. We wish to derive an equation which could serve as an approximation to the master equation 
^ in the diffusion limit. We take the Fourier transform of the Eq.Q and insert Q{k) in the above form, that yields 

^^M ^ -K'^\k\'^THx)p{x, t)]. (13) 

The Ea. ljl^ll can be formally inverted by expressing the invert Fourier transform by a fractional Weyl-Riesz operator: 
J^~^{—\k\^^) ~ g?u [^. The resulting equation is the following 

-^r^"^ — ^A- — ■ ('') 

Our aim is to solve the Ea. H14|l for v given by Eq.JHl), where 9 > —1, with the initial condition p{x,0) = S{x). 
Obviously, the presence of the x-dependent term under the fractional derivative poses the main difficulty. 

A. The case 6 = 

The Ea. (|13(l for this case can be easily solved: 

p{k,t) = exp{-K''t\k\''). (15) 

The above expression is the characteristic function of the Levy process [Sj and its inversion produces the symmetric 
Levy distribution. To handle those distributions, it is convenient to apply Fox functions. In fact, the Levy distribution 
in its most general form can be expressed as a Fox function 31J. That formalism is useful in describing stable processes 
because it reflects scaling properties of the underlying phenomena. The definition and some properties of the Fox 
functions are presented in Appendix. 

The solution of the Eq. H14(l can then be expressed in the following way 0, |23 



//|x| ■ 



{KHy/t^ 



(1,1^), (1,1/2) 
(1,1), (1,1/2) 



(16) 



The asymptotics |a;| -^ oo results from the expansion (jA3|l : p{x,t) r^ f/|a;|^+'^ and implies that the variance, as well 
as all fractional moments of the order of /i or higher, diverge. 

B. The general case 

The Ea. H14|l . as an approximation to the master equation in the diffusion limit, has been obtained from the Eq.Q 
by dropping all terms of the order higher than fc'^ in the expansion of Q{k). Consequently, solving the Ea. H14|) . we 



can neglect those terms without introducing any additional idealization. We look for a solution in the form of the 
Fox function 



p{x, t) = NaH^-j 



(ai,Ai), (02,^2) 

ibi,Bi),{b2,B2) 



(17) 



similar to that for the case 6 — 0, where a — a{t) and N is the normalization constant. We want to determine the 
coefficients a^, Ai, hi and Bi, as well as the function a(i), by expansion of both sides of the Ea. ll3|) in consecutive 
fractional powers of k and neglect the terms which are small for k ^ Q. 

The function po = x^^p{x,t) on the right hand side of Ea. H13l) can be expressed as the Fox function due to the 
multiplication rule (|A5(I . Its Fourier transform, given by the Ea. (|A6|l . expanded according to the formula l|A3|l . takes 
the form 



,W 



pg{k, t) = 27r{h\;> + hf'k' + 0{k^)), 



(s). 



(18) 



where k' = Kk/a. Deriving explicitly consecutive terms and utilizing simple properties of the gamma function, one can 
find conditions for the coefficients. First, to get the term of order fc°, we introduced the condition 1 — oi = Ai{l — 0). 

have: 



/li"'' = identically because the gamma function in the denominator is infinite. Similarly, for the function p{x,t) we 



p{k, t) = 27r(/io + h_gk'~'^ + hik' + h2k'' 



h^k'^ 



0{k'^+^)), 



(19) 



where we imposed the condition 2 — ai = Ai(l + /i) to get the exponent /i for the term A;^. Then the above conditions 
determine the first two coefficients of the Fox function: ai and Ai. We need also h-g = 0; this requirement can be 
satisfied if the argument of one of the gamma functions in the denominator assumes an integer and non-positive value. 
The condition for that is 1 — 62 ~ -82(1 — 6) — (alternatively, we could impose a similar condition for 02 and ^2)- 
The same procedure allows us to satisfy the requirement h2 — Q and to determine B2 = 1/(2 + 0). Finally, since the 
coefficients 02 and A2 enter the above expressions in a similar way as 62 and i?2j we want to preserve this symmetry, 
present for the case 6* = 0, and assume 02 — 62 and A2 — B2- 

Then we insert the expansions Eas. (|18() and H19() into the Ea. (|13|l and separate the time-dependent term: 
{a/K)^^^^^^d/K = —A, where A is a constant which scales the time and then it is not essential. We assume 
X = l/{lJ, + 9). By taking into account the initial condition for the Ea. (|14|l . one can write down the solution for the 
function a{t) as: 



a = Kt f'+s 



(20) 



and then reduce the problem to a simple equation: 



M 



'0 1 



(21) 



where h, 



W 



-BuW 



After implementing the coefficients we have evaluated, we obtain for the term hfj, the following expression: h^ = 
7r^^(^ + 9)T{—fj,)r{bi + Bi{l + /i))cos(/i7r/2)sin(^ig7r). Unfortunately, the term h^ cannot be obtained directly 

(g\f 

from the series expansion because the undetermined expression emerges. Instead, we proceed as follows, h^ ' can be 
expressed as: /ij, — {2Tr)^ ^a~^pg{k = 0) = n^^ J^ z~^W{z)dz, where we obtained the function 



W{Z) = <2 



(.^' L1+0)' (.-'-' 2+e) 



ti+0 
(6i+Bi(l 



),5i),(l 



' 2+e 



(22) 



by applying the relation (|A5|) . The required term can now be easily evaluated as the Mellin transform J^{W){s) = 
/„°° W{z)z'-^dz: 



,W' = ^-^M{W){s = 0) = ^-^(0) = 1^-' limx(6) = -^nbi + B,{1 - 9)). 



(23) 



Inserting the expressions for /i^ and % to the Ea. (|21|) . yields the relation between bi and Bi. Then, finally, the 



solution of the Eq. (|14(l reads 



p{x, t) = NaH2 2 






(6i,Si),(l - Irg, 2x5) 



with the condition 



^^r(-.TO + B,a + .„sin(|±^. 



TT cos(/i7r/2) + r(6i +Bi{l- 6)) = 



(24) 



(25) 



which follows directly from the Ea. (|21l) . The general theory of the Fox functions implies that p{x,t) is the analytic 
function for all x ^ ii Bi > l/(/i + 0). Apart from that -very weak- condition, one parameter is arbitrary. We will 
return to this issue in Subsec. D. The normalization factor N can be determined in a simple way by using the Mellin 
transform: N = [2x(~l)]~"'^, that yields 



--I 



r(6i + Bi)r - 



M- 



sm 



2 + 1 



In the fc-space, the solution (|24|) is of the form 

p(fc, t) = l- cr^fc^' + • • • « exp(~cr^fc'') 
where 



a'^ 



H 2 + e T{h,+B,)T{-^,) 



ti^+o 



(26) 



(27) 



(28) 



and we neglected the terms of the order A:^^+^ {6 < 0) or fc^^ {6 > 0). The formula (|27|l means that within the scope 
of validity of our approximation the solution H24f) is equivalent to the Levy stable distribution. 

The Fox functions can be evaluated from series expansions. Since they are poorly convergent, one needs the series 
both for small and large values of the argument. According to Ea. ljA3|) . the expansion oi p{x, t) for small |a;| is of the 
form 



p(x,t) = x''i/'^i Vr — 



1 bi+i 



ttBi ^—' \u + 9 n + 9 Bi 

Using the property (|A4() . we obtain the series for large \x 
Na 



sm 



1-1 



1 bi+i 
2 + e Bi 



{-iyia\x\y^^'/i\. (29) 



p{x,t) 



V . . / -^ ft \ 

-^{^i + 0) ^ T[b^ + B^{l-e + ^{^i + 6))] sin ( f^*^ j (-l)^(a|xr i+^^^^+^V^!- 



(30) 



The above expression implies that the tail of the distribution has the same x-dependence as for the case 9 — Q: 
p{x,t) - iM/(M+e)/|a.|i+A' Qx\ -^ cx)). 



C. Diffusion 



A characteristic feature of the Levy distributions is the divergence of moments. In particular, the mean squared 
displacement is infinite for any time and then the transport phenomena require a modified formalism for the diffusion. 
One possibility is to substitute the variance by a moment of the order 5 < ji. 

In order to evaluate the moments of the distribution H24|) . we utilize properties of the Mellin transforms. A simple 
calculation yields: 



(|a;n = 2N x^p{x, t)dx = 2Na-''x{-6 - 1) 
Jo 



2NK-'' 



T{b,+B,il + S))T 



S 



n + e 



i^Vv(M+e). 



(31) 



Applying the above expression, one can compare individual cases in respect to the speed of transport. However, as 
long as the parameter S is arbitrary, such formalism seems to be incomplete. Can it be fixed in some way? Clearly, 
the value S = fi is distinguished. Since that moment is divergent, we consider S = fi — e^ where < e <^ and then 
the case 6* = is excluded. In the limit e ^ 0, the gamma function can be expanded and Ea. (|31|l takes the form 



r-) - ^^(^ + (^nh + i?i(l + M)) Bin (f±^vr) t-/(-+«). 



{\x 



Let us now define the fractional diffusion coefficient V^f^^t): 



j)(p) = 



r(I + fi) t e-,0+ 



\ime{\xr^) (t^cx)). 



(32) 



(33) 



where f < /i < 2. According to 1(22), the limit is finite and P^^' -- ^-^/(M+e)^ 

The interpretation of the above result is straightforward. If 6* < the coefficient V^'^^ rises with time to infinity and 
we have the"superdiffusion". Conversely, for 6* > there is the " subdiffusion" . Therefore we have obtained formally 
the same result as for FPE, Ea. (|ll|l . except the variance has been substituted by the fractional moment. The common 
conclusion, drawn both from the Gaussian and the fractional case, is that the sign of 9 decides which kind of diffusion 
the system reveals. 

D. Numerical examples 

In this Subsection we evaluate the probability distributions for specific cases. They are compared with the solutions 
of the master equation Q, obtained by the Monte Carlo simulations of stochastic trajectories of the jumping process. 
The Levy-distributed jump-size density has been generated by using the algorithm from Ref . |33| . 

It follows from the expansion H29|l that p{x,t) ~ |a;|''i/'^i for small |x|. Therefore, the ratio of the coefficients bi 
and Bi determines the shape of the distribution p{x,t) there. Since our approximation is suited for large |a;|, this 
ratio remains undetermined. The analysis of the master equation indicates that its solutions exhibit the power law 
dependence of the form |a;|^ for small |a;|. We utilize this property and assume the relation bi — 6Bi in the present 
Subsection; Bi follows from the numerical solving of the Ea. (|25|l . The Fox functions have been evaluated from the 
series H29|l and H3()|l with sufficient precision for all |x| G (0, oo). 







FIG. 2: The solutions of Eq.lI^J for 6 = —0.2 (dots) and — 1 (dashed line). The corresponding solutions 12411 are marked by 
solid lines. 



The solutions of the fractional equation (|24|) . for a negative and a positive value of 9, are presented in Fig. 2. They 
correspond to the superdiffusion and subdiffusion, respectively, and display very different shapes. All the distributions 



8 

with 9 > rise at small \x\ and display a maximum, whereas those for 6 < fall monotonically. The comparison 
with the solutions of the master equation is also presented in Fig. 2. The curves which represent both equations have 
similar shapes and their tails coincide. 

Fig. 3 presents the comparison of the Fourier transforms for the solution of both equations; the same cases as in 
Fig. 2 are shown: 9 = —0.2 and 9 ~ I. For that purpose, the master equation has been solved for all x up to very 
large values, then the numerical integration with the cosine function has been performed. In the case of the master 
equation solutions, the figure reveals substantial deviations from the straight lines - which represent the shape of the 
Levy distribution - for large k: for 9 < that distribution stabilizes with fc, whereas it falls rapidly to zero for 9 > 0. 
At small k both solutions coincide with those of the fractional equation. 







0.01 



0.001 



0.0001 : 



FIG. 3: Logarithm of the Fourier transform from the solution of Eq.(|5J for 9 — —0.2 (dots) and 9 — 1 (dashed line). The 
same quantity for the solution H24^ which corresponds to both cases, calculated according to Eg. 127^ . is marked by solid lines; 
^ = k/a. 



V. DISCUSSION 



The jumping process we have discussed in this paper is Markovian because the waiting time probability distribution 
has the exponential form. However, since the jumping frequency depends on the process value, the system possesses 
some properties which are usually attributed to non-Markovian processes. In particular, we have demonstrated the 
presence of the anomalous diffusion. 

If the step-size is small, the jumping process can be regarded as continuous and described in terms of the FPE. 
That approximation has been accomplished by applying the Kramers-Moyal expansion, on the assumption that all 
moments of the step-size distribution are finite. It has been demonstrated -by solving the FPE exactly for the 
power-law frequency i^{x)- that both normal and anomalous diffusion emerges. 

On the other hand, we have considered the diffusion limit of the master equation for the step-size distribution which 
is stable and has divergent moments. The resulting equation (I14() is fractional and contains a variable coefficient, 
in contrast to usually studied equations which govern Levy processes. We have demonstrated that in the diffusion 
limit the Ea. H14(l is satisfied by the Fox function if ly depends algebraically on the position. The coefficients of the 
Fox function have been derived by inserting it to the equation and by comparing the terms. However, one parameter 
must remain undetermined if only the diffusion limit is considered because it is responsible for the behaviour of the 
solution at small |x|. 

The solution, since it is expressed as the Fox function, depends on time in the scaling form. Due to that property, 
simple conclusions about the transport can be drawn. The fractional moments are given by complicated expressions 
but the time-dependence factorizes and it obeys a simple power law. We have defined the fractional diffusion coefficient 



PCa") which allows us to establish a correspondence with the standard description in terms of the FPE. Both approaches 
predict the diffusion coefficient, either T) or 'D^'^\ in the form ~ ^-^/(a'+s)^ [q tjie subdiffusion {9 > 0) and the 
superdiffusion {6 < 0). We would like to emphasize that the above conclusions are independent of a specific choice of 
the free parameter in the solution. 

An independent numerical analysis of the master equation (0) provides additional information about the jumping 
process. We can learn that it is not the Levy stable process: substantial deviations from the Levy distribution at 
large k are clearly visible. They become meaningless in the diffusion limit. 

The comparison of the solutions of both diffusion equations we have discussed in this paper with the results of the 
numerical analysis of the master equation shows a reasonable agreement not only in the asymptotic limit of large \x\. 
Since the expressions (O and (|24|l are relatively simple, compared to the Eq.©, both diffusion equations could serve 
as convenient approximations to the master equation |(2Jl. 



APPENDIX 



The Fox function |34L \33L l3q| (for a review of the most important properties see [iJi 1^ ) is defined as an inverse 
Mellin transform in the following way: 



Trmn 

pq 



[ttp, Ap) 



Tjmn 



{ai, Ai),{a2, A2), . . . ,{ap, Ap) 
ibi,B,),{b2,B2),...,ib„B,) 



2TTi 



x{s)z'^ds, 



where 



Xis) 



UTnbj-B,s)Uim-a,+A,s) 

Ul+1 r(i - b, + B,s) nr,+i r(a, - A,s) ■ 



(Al) 



(A2) 



Coefficients Ai and Bi are positive. The contour L is a straight line parallel to the imaginary axis which separates 
the poles of both gamma functions in Ea. (|A2(l . If those poles are simple, the Fox function can be expressed in the 
form of the following series 



Trmn 

pq 



[ap, Ap) 
iKBq) 



m 00 



y y il7=l,,^. r(&, - i?,^) n;Li r(l - a, + A, ^) (^l).,(b.+.)/B. 



A similar expansion can be obtained for z — * 00 by using the property: 



Trmn 

pq 



[Up, Ap) 
ibq,Bg) 

Another useful property is the multiplication rule: 

(flp, Ap) 
iKB,) 

The Fourier cosine transform of the Fox function yields: 

(flp, Ap) 



a Trmn 



Trmn 
-iipq 


1 

Z 


{l-ap,Ap) _ 




r 


{ap + aAp,Ap) 


^Bpq 


Z 


{bq+aBq,Bq) 



Tj-rnn 



ibq,B 



qj -^qJ 



cos{kx)dx = -r^^+i,p+2 



k 



il^bq,Bq), (1,1/2) 

(l,l),(l-ap,Ap), (1,1/2) 



(A3) 



(A4) 



(A5) 



(A6) 
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